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The Drell-Yan differential cross section is measured in pp collisions at y's = 7 TeV, from a data sample collected 
with the CMS detector at the LHC, corresponding to an integrated luminosity of 36 pb~^. The cross section 
measurement, normalized to the measured cross section in the Z region, is reported for both the dimuon and 
I dielectron channels in the dilepton invariant mass range 15-600 GeV. The normalized cross section values are 

I quoted both in the full phase space and within the detector acceptance. The effect of final state radiation is 

«^ also identified. The results are found to agree with theoretical predictions. 

(N 

<U 1. Introduction 

The production of lepton pairs in hadron-hadron collisions via the Drell-Yan (DY) process is described in the 
standard model (SM) by the s-channel exchange of 7*/Z. Theoretical calculations of the differential cross section 
da / dM {££) , where M{££) is the dilepton invariant mass, are well established up to next-to-next-to-leading order 
^ (NNLO) [TH5|. Therefore, comparisons between calculations and precise experimental measurements provide 
O stringent tests of perturbative quantum chromodynamics (QCD) and significant constraints on the evaluation of 
^ ^ the parton distribution functions (PDFs). Furthermore, the production of DY lepton pairs constitutes a major 
source of background for top quark pair and diboson measurements, as well as for searches for new physics, 
such as production of high mass dilepton resonances. 

Here we present a measurement of the differential DY cross section in proton-proton collisions at ^/s = 
I 7 TeV, based on dimuon and dielectron data samples collected in 2010 by the Compact Muon Solenoid (CMS) 
^ experiment at the Large Hadron Collider (LHC), corresponding to an integrated luminosity of 35.9 ± 1.4 pb~^. 
l/^ The results are given for the dilepton invariant mass range 15 < M{££) < 600 GeV. In the analysis presented, 
the cross sections are calculated as 

I where is the unfolded background-subtracted yield, corrected for detector resolution. The values of the 

't^ acceptance A and the efficiency e are estimated from simulation, while p is a factor that accounts for differences 
• ^ in the detection efficiency between data and simulation. Knowledge of the integrated luminosity C is not 
required for the measurements described in this paper, since the cross sections are normalized to the Z region 
^ (60 < M{££) < 120 GeV). The measurements are described in detail in Ref. g]. 



2. The CMS Detector 

A detailed description of the CMS detector and its performance can be found in Ref. 53 . The central feature 
of the CMS apparatus is a superconducting solenoid 13 m in length and 6 m in diameter, which provides 
an axial magnetic field of 3.8 T. Within the field volume are the silicon pixel and strip tracker, the crystal 
electromagnetic calorimeter (ECAL), and the brass/scintillator hadron calorimeter (HCAL). Charged particle 
trajectories are measured by the tracker, covering the full azimuthal angle and pseudorapidity interval \r]\ < 2.5, 
where the pseudorapidity is defined as 77 = — lntan(6'/2), with 6 being the polar angle of the trajectory of the 
particle with respect to the counterclockwise beam direction. Muons are measured in the pseudorapidity range 
hi < 2.4. 



3. Event Selection 

The analysis presented in this paper is based on dilepton data samples selected by inclusive single-lepton 
triggers with px thresholds ranging between 9 and 15 GeV for muons and between 15 and 17 GeV for electrons, 
depending on the beam conditions. 
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Figure 1: The observed dimuon (left) and dielectron (right) invariant mass spectra. No corrections are applied to the 
distributions. The points with error bars represent the data, while the various contributions from simulated events are 
shown as stacked histograms. By "EWK" we denote W ^ £v and diboson production. The "QCD" contribution results 
from processes associated with QCD and could be genuine or misidentified leptons. The lower panels show the ratios 
between the measured and the simulated distributions including the statistical uncertainties from both. 



Muons are required to pass the standard CMS muon identification and quality criteria, based on the number 
of hits found in the tracker, the response of the muon chambers, and a set of matching criteria between the 
muon track parameters as determined by the inner tracker section of the detector and as measured in the muon 
chambers [51 [7] . Cosmic-ray muons are further suppressed by requirements on the proximity to the interaction 
region and the opening angle between the two reconstructed muons. The two muons are fitted to a common 
vertex reducing the contamination from QCD related events. Finally, an isolation requirement is imposed on 
each muon, based on the relative transverse energy deposited in the HCAL and transverse momenta of tracks 
around each of the muons. 

Electron reconstruction starts from clusters of energy deposited in the ECAL, and associates with them hits 
in the CMS tracker [8]. Energy-scale corrections are applied to individual electrons as described in Ref. [9]. 
The electron candidate is required to be consistent with a particle originating from the primary vertex in 
the event. Electron identification criteria based on shower shape and track-cluster matching are applied to 
the reconstructed candidates. Electrons originating from photon conversions are rejected by eliminating those 
electrons for which a partner track consistent with a conversion hypothesis is found, and requiring no missing 
hits in the pixel detector, as discussed in Rcf. |10| . Isolation requirements are imposed on each electron, based 
on the relative transverse energy deposited in the ECAL and the HCAL and transverse momenta of tracks 
around each of the muons. 

Muons must be reconstructed with \r]\ < 2.4 and px > 7 GeV but at least one of them is required to be with 
\r]\ < 2.1 and px > 16 GeV in order to ensure effective triggering of the event. Electrons must be reconstructed in 
the ECAL barrel with \ri\ < 1.44 or in the ECAL endcaps with 1.57 < \r]\ < 2.5. The leading electron is required 
to have Ex > 20 GeV providing it triggered the event, while the second electron must have Ex > 10 GeV. 

Event samples for simulation studies of electroweak processes involving W and Z production are produced with 
the NLO MC generator POWHEG [TTHTg] interfaced with the PYTHIA (v. 6.422) [H] parton-shower event 
generator, using the CTIO [1.5. parametrization of the PDFs. PYTHIA is also used for the FSR simulation. 
The QCD multijet background is generated with PYTHIA, and background from top quark pairs is simulated 
using MadGraph (v. 4.4.12) [TS] and PYTHIA at leading order using the CTEQ6L PDF set [H] for both 
samples. Generated events are processed through the full GEANT4 |TH] detector simulation, trigger emulation, 
and event reconstruction chain. 

The observed invariant mass distributions, in the dimuon and dielectron channels, are shown in Fig. [T] 
Thirteen mass bins with unequal widths are used to cover the observable dilepton mass spectrum. 
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4. Backgrounds 

The main backgrounds at high dilepton invariant masses are caused by top quark pairs and diboson pro- 
duction, while at invariant masses below the Z peak, DY production of r+r^ pairs becomes the dominant 
background. At low dimuon invariant masses, most background events are QCD multijet events that originate 
from both heavy flavor and light quark production. The expected shapes and relative yields of these several 
dilepton sources can be seen in Fig. [T] They are in agreement with the estimates extracted from data and the 
significant background sources are estimated from datajl]. 



5. Detector Resolution Effects and Unfolding 

The effects of the detector resolution on the observed dilepton spectra are corrected through an unfolding 
procedure. The original invariant mass spectrum is related to the observed one (in the limit of no background) 
by 

^obs,i = ^ Tik Ntruc,k, (2) 
k 

where A''^ is the event count in a given invariant mass bin i. The element Tik of the "response matrix" T is 
the probability that an event with an original invariant mass in the bin k is reconstructed with an invariant 
mass in the bin i. The original invariant mass spectrum is obtained by inverting the response matrix and 
calculating [50] 

Nu,k = -^truo,fe = ^ {T ^)ki Nohs,i- (3) 
i 

This procedure is sufficient in the analysis reported in this paper because the response matrix is nonsingular 
and nearly diagonal. Two extra dilepton invariant mass bins are included in the unfolding procedure, to account 
for events observed with M{U) < 15 GeV or M{ee) > 600 GeV. 

The response matrix is calculated using the simulated sample of DY events, defining the "true mass" as the 
"generator level" dilepton invariant mass, after FSR. Only the selected events in the sample are used to calculate 
the response matrix. The loss of events caused by reconstruction inefficiencies or limited acceptance is factored 
out from the unfolding procedure and taken into account by means of efficiency and acceptance factors in a 
subsequent step. The response matrices in both lepton channels are invertible and are represented on Fig. |2] 




generated M(n|i) (post-FSR) [GeV] generated M(ee) (post-FSR) [GeV] 

Figure 2: The response matrices for the muon (left) and electron (right) channels from simulation. 
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6. Acceptance and Efficiency 

The geometrical and kinematic acceptance A is defined, using the simulated leptons after the FSR simulation, 
as A = -/VaccZ-^gonj whcrc A'gon is the number of generated events and iVacc is the corresponding number of events 
passing the standard px and r] lepton requirements, in each dilepton invariant mass bin. 

The efficiency e is the fraction of events within the acceptance that pass the full selection, so that 



A-e. 



(4) 



where is the number of events surviving the reconstruction, selection, and identification requirements. The 
values of the product of acceptance and efficiency are obtained from simulation. A separate correction factor is 
determined from data and applied to the product, following the procedure used in the inclusive W and Z cross 
section measurements in CMS ^10]. This factor, the efficiency correction, describes the difference between data 
and simulation in the efficiency to observe single leptons or dileptons. 

The POWHEG simulation combines the next-to-leading-order (NLO) calculations with a parton showering 
which is insufficient to model fully the low invariant mass region of the dilepton spectra. The two high-px leptons 
required in the analysis must form a small angle at low mass and therefore the dilepton system gets significantly 
boosted, something to be compensated by hard gluon radiation in the transverse plane. This means that these 
low-mass events are of the type "7* -I- hard jet" at first order, and therefore the next order of correction (NNLO) 
becomes essential for a reliable estimate of acceptance corrections. To account for this, a correction is applied, 
determined from the ratio between the differential cross sections calculated at NNLO with FEWZ [2Tj and at 
NLO with POWHEG, both before FSR. These correction weights, obtained in bins of dilepton rapidity, px, 
and invariant mass, are applied on an event-by-event basis. This procedure changes the acceptance in the lowest 
invariant mass bin significantly (by about 50%), but has a small effect, not exceeding 3%, on the rest of the 
bins. 

Figure [3] shows the variables A, e, and A ■ e as functions of M{££) for dimuons (left) and dielectrons (right). 

The FSR correction factors for a given invariant mass range are obtained from simulation by dividing the 
cross sections after FSR by the corresponding quantities before FSR. The factors are channel and invariant mass 
dependent and are up to few hundred of percent just below the Z peak but 10-20 % or less elsewhere. They are 
applied on (corrected) data as an additional step. The factors obtained within the detector acceptance and in 
the full phase space are applied to the corresponding measurements. 
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Figure 3: DY acceptance (blue, filled circles), efficiency (red, open triangles), and their product (black, open squares) 
per invariant mass bin, for the fJ.'^fi' (left) and e''"e~ (right) channels. 

The total dimuon event selection efficiency is factorized as 



£(event) = e(^i) • £(^^2) ■ e[MMl(Mi)&(M2)] • £(event, trig|^^). 



(5) 



where e(/i) is the single muon selection efficiency; £[/Lt/Lt|(/Lii)&(/Lt2)] is the dimuon selection efficiency, which 
includes the requirement that the two muon tracks be consistent with originating from a common vertex and 
that they satisfy the angular criteria; and e(event, trigj/i/j,) is the efficiency of triggering an event. The single 
muon efficiency is factorized as 
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= e(track| accepted) • e(reco + id|track) • e(iso|reco + id), (6) 

where e(track|accepted) is the offline track reconstruction efhciency in the tracker detector; e(reco + id | track) 
is the muon reconstruction and identification efficiency; and £(iso|reco + id) is the muon isolation efficiency. The 
trigger efhciency e(event, trig|/x/i) is given by 



£(event,trig|^/i) = e(/^i, trig|^i) + e(A*2, trig|//2) - £(//i, trig|^i) • e(^2, trigj^a), (7) 
where e(/i, trig|/z) is the efhciency of an ofhine selected muon to fire the trigger. 

The track reconstruction efficiency is very high (99.5%). The angular criterion is nearly 100% efhcient for 
signal DY events, and the vertex probability requirement is more than 98% efficient and has a negligible (< 0.3%) 
dependence on M{££). 

The muon reconstruction and identification efficiency is estimated using clean samples of muon pairs in the Z 
peak (tag and probe, T&P, method [lO]). To determine the isolation efficiency, the Lepton Kinematic Template 
Cones (LKTC) method IS] is applied. 

To describe the observed efficiency variations between data and simulation, efficiency correction factors are 
obtained in bins of px and rj as the ratio of the efficiencies measured with data and with the simulated events: 

Pcff PT , V) = -, r ■ (8) 

The corrections to the efficiencies in simulation are implemented by reweighting simulated events, with weights 
computed a,s W = p\°'^°^^'^ P2°'^°^^'^ pf° P2° P*""^ where pf is estimated per muon i — 1,2 from Eq. using 



^ - "reco + ^d","z.o" efiiciencies defined in Eq. ^ and = (4aL,i + 4ata,2 " 4afa.i4aL,2)/(^m,i + 
^s^m 2 ~ ^sim i^m 2)- PT < 16 GcV Or \ri\ > 2.1 for a given muon, its trigger efficiency is set to zero. 

The total event efficiency in the dielectron channel analysis is defined as the product of the two single electron 
efficiencies, which incorporate three factors: 1) the efficiency Ereco to reconstruct an electron candidate from 
an energy deposit in the ECAL; 2) the efficiency Sid for that candidate to pass the selection criteria, including 
identification, isolation, and conversion rejection; 3) the efficiency etrig for the leading electron to pass the trigger 
requirements. Each of these efficiencies is obtained from simulation and corrected by poff (pT) ??), as for the muon 
channel (Eq. Q). The T&P method is used for all efficiency components. The event efficiency correction and 
its uncertainty are derived as for the muon channel by reweighting simulated events. 

The overall efficiency correction is up to 10 percent with largest values at lowest invariant masses. 



7. Systematic Uncertainties 

The acceptance-related uncertainty resulting from the knowledge of the PDFs is estimated using PYTHIA 
with the CTEQ6.1 PDF set by a reweighting technique [22]. It is approximately the same for the dimuon and 
dielectron channels and is between one and three percent in the whole invariant mass range under study. 

The uncertainty of the acceptance is estimated, for each dilepton invariant mass bin, using FEWZ, at NLO 
and NNLO accuracy in perturbative QCD. Variations of the factorization and renormalization scales lead to a 
systematic uncertainty smaller than 1% (at NNLO) for most of the invariant mass range used in the analysis 
presented here. 

Since the POWHEG MC (NLO) simulation, modified to match the FEWZ (NNLO) calculations, is used 
to calculate the acceptance corrections used in the analysis, an additional (model-dependent) systematic un- 
certainty on the acceptance calculation is determined from the observed differences in acceptances based on 
FEWZ spectra and POWHEG distributions matched to FEWZ. This systematic uncertainty reaches up to 
10% (at lowest invariant masses) in the dilepton invariant mass range considered in the analysis and is included 
in the comparison between the measurements and the theoretical expectations. 

The dominant systematic uncertainty on the cross section measurement in the dimuon channel is the uncer- 
tainty on the background estimation. It reaches four percent at very low invariant masses and is also substantial 
at high masses where the measurement is statistically limited. The next most important uncertainties are re- 
lated to the muon efficiency and to the muon momentum scale and resolution, directly related to the unfolding 
procedure. Dimuons from the Z region and the invariant mass shape distribution are used to assess these uncer- 
tainties. They are of the order of one percent except in the proximity to the Z peak (2-4 %) and also at very high 
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masses where statistics is poor ( 2%). Comparisons of the FSR spectrum between data and simulation show 
that potential differences could not affect the cross section measurements by more that 2% and are typically 
confined within 0.5 %. Remaining systematic effects are within f % and/or not exeeding the statistical precision 
of the measurement. 

In the electron channel, the leading systematic uncertainty is associated with the couple of percent uncertainty 
on the energy scale correction of individual electrons. Its effect on the cross section measurement is very 
significant and reaches tens of percent at very low invariant masses and around the Z peak, it is several percent 
elsewhere. The second leading uncertainty for electrons is caused by the uncertainty on the efficiency scale 
factors which is statistically limited. It is close to 10% in the low invariant mass region and few percent in the 
rest of the spectrum. The dielectron background uncertainties follow the pattern described for muons but are 
slightly higher at low invariant masses reaching 4-6 %. The resolution effects through the unfolding procedure 
are smaller than the leading source for each of the invariant mass bins and are typically few percent. Because of 
significantly higher systematic uncertainty for all mass bins for the electron channel than for the muon channel, 
the FSR related contribution to the electron channel systematic uncertainty is neglected. 



8. Results 

The DY cross section per invariant mass bin i, cr^, is calculated according to Eq. ([T]). 

In order to provide a measurement independent of the luminosity uncertainty and to reduce many systematic 
uncertainties, the ai is normalized to the cross section in the Z region, au, defined as the DY cross section in 
the invariant mass region 60 < M{££) < 120 GeV. The result of the analysis is presented as the ratio 

ni _ ^u,i I -^n.norm /„v 

-«post-FSR - A I 4 - > W 



^norm ^norm Pnorin 



where A^u,i is the number of events after the unfolding procedure, and the acceptances Ai, the efficiencies e^, 
and the corrections estimated from data, pi, were defined earlier; A^u^norm, ^norm, ^norm, and Pnorm refer to the 
Z region. For both lepton channels, the cross sections in the Z region measured in this analysis are in excellent 
agreement with the previous CMS measurement [TD]. 

In order to allow a more direct and precise comparison with theory predictions, the shape measured before the 
acceptance correction is also reported, thus eliminating PDF and theory uncertainties from the experimental 
results: 

r>i _ -^u,i / -^u.norm 

-ftdet, post-FSR — ^ „ / p „ ■ y^^l 

Pi '^norin pnorm 

The shapes corresponding to the DY process after FSR, -Rpost-FSR and -Rdot,post-FSRi are modified by the FSR 
correction factors, to obtain the shapes before FSR, R and i?dcti respectively. The shapes integrated in the 
normalization region are equal to one by construction. 

The results are presented in Tables [l] and [ll] respectively, for the dimuon and dielectron channels. The two 
shape measurements, shown in the last column of the tables, are in good agreement for 11 out of 13 invariant 
mass bins and remain statistically consistent (although marginally) for the remaining two bins, 40-50 GeV and 
120-150 GeV. 

The results are also normalized to the invariant mass bin widths, AM^, defining 

Assuming lepton universality, the dimuon and dielectron results for are combined in a weighted average, 
using as weights the inverse of the respective squared total uncertainties, where the statistical and systematic 
uncertainties are added in quadrature. 

Figure |4] compares the measured (combined) results for the shape r with the prediction from the FEWZ 
NNLO calculations, performed with the MSTW2008 PDF set [23]. The measurements are very well reproduced 
by the theoretical calculations. 

More details on the measurements can be found in Ref . [3] . 
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Table I: Results for the DY spectrum normalized to the Z region in the dimuon channel. The statistical and systematic 
uncertainties are summed in quadrature. i?post-FSR and J?dot,post-FSR are calculated using Eqs. (|9| and (10 1, respectively. 
The -Rdct and R are calculated using the FSR corrections. 



Invariant mass bin (GeV) 


-Rdot.post-FSR (10 ^) 


i?dct (10-^) 


J?post-FSR (10 ^) 


R (10 




15-20 


18 


± 


2 


19 ± 


2 


772 


± 


67 


780 ± 


69 


20-30 


58 


± 


3 


58 ± 


3 


528 


± 


33 


533 ± 


34 


30-40 


67 


± 


3 


67 ± 


3 


147 


± 


8 


147 ± 


8 


40-50 


44 


± 


2 


41 ± 


2 


66 


± 


4 


62 ± 


4 


50-60 


30 


± 


2 


23 ± 


2 


37 


± 


3 


30 ± 


2 


60-76 


51 


± 


2 


28 ± 


1 


55 


± 


3 


32 ± 


2 


76-86 


97 


± 


4 


56 ± 


3 


98 


± 


5 


58 ± 


3 


86-96 


803 


± 


14 


861 ± 


15 


799 


± 


23 


857 ± 


26 


96-106 


38 


± 


3 


43 ± 


3 


37 


± 


3 


41 ± 


3 


106-120 


12 


± 


1 


12 ± 


1 


11 


± 


1 


12 ± 


1 


120-150 


9.2 


± 


0.9 


9.7 ± 


1.0 


8.4 


± 


0.8 


8.8 ± 


0.9 


150-200 


3.1 


± 


0.6 


3.2 ± 


0.7 


2.6 


± 


0.5 


2.7 ± 


0.6 


200-600 


1.8 


± 


0.4 


1.9 ± 


0.5 


1.4 


± 


0.3 


1.5 ± 


0.4 



Table 11: Results for the DY spectrum normalized to the Z region in the dielectron channel. The statistic al a nd systematic 
uncertainties are summed in quadrature. i?post-FSR and -Rdot,post-FSR are calculated using Eqs. ||9| and (10 1, respectively. 
The i?dot and R are calculated using the FSR corrections. 



Invariant mass bin (GeV) 


-Rdct.post-FSR (10 ^) 


i?dot (10-=*) 


-Rpost-FSR. (10 ^) 


R (10^=*) 


15-20 


6 


± 


3 


6 ± 3 


487 


± 


230 


508 


± 238 


20-30 


13 


± 


2 


13 ± 2 


536 


± 


96 


559 


± 97 


30-40 


24 


± 


4 


22 ± 4 


129 


± 


22 


131 


± 21 


40-50 


28 


± 


4 


24 ± 4 


52 


± 


8 


47 


± 7 


50-60 


30 


± 


5 


19 ± 3 


39 


± 


6 


27 


± 4 


60-76 


78 


± 


12 


30 ± 4 


84 


± 


13 


36 


± 5 


76-86 


144 


± 


60 


61 ± 25 


147 


± 


60 


64 


± 26 


86-96 


722 


± 


62 


839 ± 60 


715 


± 


62 


834 


± 60 


96-106 


44 


± 


21 


55 ± 26 


43 


± 


20 


53 


± 25 


106-120 


13 


± 


3 


15 ± 3 


12 


± 


2 


14 


± 3 


120-150 


5.4 


± 


1.2 


6.0 ± 1.3 


4.8 


± 


1.1 


5.4 


± 1.2 


150-200 


2.5 


± 


0.8 


2.8 ± 0.8 


2.1 


± 


0.6 


2.3 


± 0.7 


200-600 


2.1 


± 


0.6 


2.4 ± 0.7 


1.5 


± 


0.5 


1.7 


± 0.5 



9. Conclusions 

The Drell-Yan differential cross section normalized to the cross section in the Z region has been measured in 
pp collisions at ^/s = 7 TeV, in the dimuon and dielectron channels in the invariant mass range 15 < M{i£) < 
600 GeV. The measurement is based on event samples collected by the CMS experiment, corresponding to 
an integrated luminosity of 35.9 ± 1.4 pb~^. Results are presented both inside the detector acceptance and in 
the full phase space, and the effect of final state QED radiation on the results is reported as well. A correct 
description of the measurements requires modeling to NNLO for dilepton invariant masses below about 30 GeV. 
The measurements are in good agreement with the NNLO theoretical predictions, as computed with FEWZ. 
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Figure 4: DY invariant mass spectrum, normalized to the Z resonance region, r — (l/au)da/dM(U), as measured 
and as predicted by NNLO calculations, for the full phase space[4]. The vertical error bar indicates the experimental 
(statistical and systematic) uncertainties summed in quadrature with the theory uncertainty resulting from the model- 
dependent kinematic distributions inside each bin. The horizontal bars indicate bin sizes and the data points inside are 
placed according to Ref. [24j . The width of the theory curve represents theoretical uncertainties which do not exceed 
few percent. 
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